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Abstract. This paper describes a general approach to the compart mental modeling of nuclear 
data based on spectral analysis and statistical optimization. We utilize the renal physiology as test 
case and validate the method against both synthetic data and real measurements acquired during 
two micro-PET experiments with murine models. 

1. Introduction. Positron Emission Tomography (PET) [15] is an imaging tech- 
nique capable of detecting picomolar quantities of a labeled tracer with temporal res- 
olution of the order of seconds. FDG-PET \W\ UM is a PET modality in which 
[ 18 F]fluoro-2-deoxy-D-glucose(FDG) is used as a tracer to evaluate glucose metabolism 
and to detect diseases in many different organs. From a mathematical viewpoint, PET 
(and, specifically, FDG-PET) experiments involve two kinds of problems. First, sig- 
nal processing techniques are applied to reconstruct the time dependence of location 
and concentration of tracer from the measured radioactivity. Second, these dynamic 
PET data can be processed to estimate physiological parameters that describe the 
functional behavior of the inspected tissues. This paper focuses on this second as- 
pect and, specifically, examines FDG-PET data from kidneys and bladder, corrected 
for radiactive decay, to recover information on renal physiology. The data processed 
within the paper have been acquired by means of a PET device for small animals 
(mice), but the mathematical procedure for the treatment of data is essentially the 
same as for PET data of humans. 

The analysis of tracer kinetics in tissues whose activities have been measured with 
FDG-PET is typically based on compartmental models [l7l|6j[T8]. This conceptual 
framework identifies different compartments in the physiological system of interest, 
each one characterized by a specific (and homogeneous) functional role (therefore a 
compartment is not necessarily represented by a specific organ or anatomical district). 
As any other tracer, FDG is injected into the system with a concentration mathemat- 
ically modeled by the so-called Time Activity Curve (TAC). The time dependent 
concentrations of tracer in each compartment constitute the state variables that can 
be determined from PET data. The time evolution of the state variables (the kinetics 
of the system) is modelled by a linear system of ordinary differential equations for the 
concentrations, expressing the conservation of tracer during flow between compart- 
ments. The (constant) micro-parameters describing the input/output rates of tracer 
for each compartment are called exchange coefficients or rate constants, represent the 
physiological parameters describing the system's metabolism and are the unknowns 
to be estimated. 
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kidneys 



Fig. 1.1. The compartmental model of renal physiology adopted in this paper. 

The particular case of renal physiology is described in Figure where: 

• kidneys include two compartments, tissue (parenchyma) and pre-urine, while 
a third compartment, the urine, represents the bladder, where tracer is accu- 
mulated; 

• the TAC for blood is determined from Regions of Interest (ROIs) of the left 
ventricle drawn on FDG-PET maps at different time steps; 

• estimates of the tracer concentration for the bladder and for the two com- 
partment system made of parenchyma and pre-urine are obtained by means 
of ROIs including bladder and kidneys, respectively; and, finally, 

• the set of unknowns is made of six exchange coefficients, also referred to as 
micro parameters (other exchange coefficients, for example the ones describing 
the circulation for the bladder back to the kidneys, are set to zero for well- 
established physiological reasons). 

Simpler models for tracer circulation inside kidneys can be found in the literature, 
involving a smaller number of exchange coefficients. While the reduction in number 
allows for an estimate of the unknown exchange coefficients by rather simple math- 
ematical procedures (such as graphical methods, [13 J, this implies a notably less 
refined description of renal physiology. 

The present paper introduces a novel method for the reduction of compartmen- 
tal models like this renal one, whose main advantages are a notable generality and 
limited computational requirements. This method is based on spectral analysis and 
statistical optimization and includes two steps. In the first step, we solve the direct 
problem of determining the explicit formal expressions of the tracer concentration in 
each compartment, under the assumption that the exchange coefficients are given. 
According to the properties of these coefficients, solutions are classified in four dif- 
ferent families. Consistently with a spectral analysis approach, the concentrations 
are expressed as linear combinations of the TAC convolved with specific exponential 
functions. In the second step, the inverse problem of determining the six unknown 
exchange coefficients is addressed by means of an 'ad hoc' implementation of Ant 
Colony Optimization [3j [20] , which is a statistical optimization algorithm inspired by 
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evolutionary strategies. In this application the functional to minimize measures the 
discrepancy between the experimental concentrations and the analytical forms pro- 
vided by the solution of the forward problem and also accounts for the blood fraction 
with which the vascular system of the kidneys is fed up. 

The paper content is organized as follows. In Section 2 we describe the Cauchy 
problem modeling the direct problem. Section 3 determines the exchange coefficients 
by means of statistical optimization. Section 4 shows some applications to synthetic 
and real PET measurements. Our conclusions are offered in Section 5. 

2. The direct problem. The state variables of the three-compartment model 
adopted in this paper are the tracer concentrations in the tissue (Ct), in the pre-urine 
(C p ) and in the urine contained in the bladder (C u ). Moreover, the kinetic process 
in the system is initialized by the TAC C&, representing the tracer concentration in 
blood. The usual conditions for compartmental analysis hold, e.g., tracer is uniformly 
distributed in each compartment at each instant [18] , diffusive effects are neglected and 
physiological processes are in a steady state. The six constant exchange coefficients 
(flux rate parameters) between compartments in contact are denoted as /c a 5, where 
the suffixes a and b denote the target and source compartment, respectively. For 
example, k up is the rate constant of tracer carried "to" the bladder pool u "from" the 
pre-urine pool p. We assume k a b > for all cases. 

Conservation of tracer exchanged between compartments leads to the following 
system of linear ordinary differential equations with constant coefficients: 



-(fat + k pt ) C t + k tp C p + k tb Cb, 
kpt Ct — (hp + k U p) Cp + kpb Cb, 



(2.1) 



kup Cpj 



with initial conditions C t (0) = C p (0) = C u (0) = and where the dependence on time 
is implicit. The first two equations in ( |2.1[ ) can be written in the compact form 



C= AC+Kft, 
with initial condition C(0) = 0. Here, in fact, 

c t 

Cr> 



(2.2) 



(2.3) 



(ht + k pt ) k tp 




—a b 


kpt ~(hp H~ k U p) 




c —d 



and 



K := 



hb 
kpb 



The definition of A corresponds to define 

a := ht ~h kpt, c := k p t, d := ktp + k U p, b := hp 
with inverse relations 



(2.4) 



(2.5) 



(2.6) 



kpt = c, k bt = a- c, k tp = b, k up = d - b. (2.7) 
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In the first part of this section we will show how the structure of matrix A is related 
to the explicit forms of tracer concentration in the different compartments. In the 
framework of compartmental analysis this issue is crucial for the study of the inverse 
problem performed in the next section. 



Let's first assume c / and 6^0. Then the solution of (2.2)-(2.7) is given by 
C t = c 1 bE 1 + c 2 (d + A 2 ) E 2 , (2.8) 



Ei 



C p = a (a + \ 1 )E 1 + c 2 cE 2 , 



I e x ^ t - T) C b (r)dT = e xa *C b , i = l,2, 
Jo 



(2.9) 



(2.10) 



where * denotes the convolution operator, the constants c\ and C2 are denned as 

-ck tb + (d + X 2 )k pb _ (a + X^ktb-bkpb 

Cl ~ (a + \ 1 )(d + \ 2 )-bc' C2 ~ (a + A 1 )(d + A 2 )-6c' [ ' ' 

and Ai, A2 are the eigenvalues of the matrix A with explicit values 



Ai 



-(a + d) + i/(a + d) 2 - 4(ad - be) 



(2.12) 



A 2 



-(a + d) - i/(a + d) 2 - A(ad - be) 



(2.13) 



We observe that both Ai and A2 are negative. The concentration C u is evaluated by 

(2.14) 



inserting ( |2.9[ ) in the third equation of ( |2.1[ ), which leads to 
C u (t) 



I C p (a)d(i= I [a (a + Ai)^i(cr) + c 2 cE 2 (a)] da. 
Jo Jo 



We utilize the identity 



I I C(r)e w ^- r) drda= - [ C(a) e w da - - [ CMdr, (2.15) 
Jo Jo ™ Jo w Jo 

where C is any continuous function and w is a real parameter and find 

^ El + c 2 f E 2 - [ Cl + c 2 f ] f C b (r) dr. (2.16) 

1 Ai A 2 J 



C u (t) _ (o + Ai) 

^1 



Ai 



We now suppose that 6 = and c ^ in the definition of the matrix A, i.e. 



-a 
c —d 



(2.17) 



In this case we set Ai = —a and A2 = —d. This corresponds to vanishing denominators 
in (2.11). Integration of the linear system leads to 



C t = X i (Ai - A 2 ) Ei 
4 



(2.18) 



and 



C p = cxi E 1 +X2 E 2 , 



where E\ and E 2 are convolutions defined as in (2.10) and 

hb 



Xi 



Ai — A2 



X2 = k pb 



ck tb 
Ai — A2 



Moreover, integrating the third equation in ( |2.1| ) leads to 



k 



up 



- E\£ + ^ £^2,6 ' 
Ai A2 



Ai Ai 



C 6 (r) dr. 



We then consider the case b ^ 0, c = 0, i.e. 

—a 6 
-d 



(2.19) 



(2.20) 



(2.21) 



(2.22) 



Again we have Ai = —a and A2 = — d. The solution of the Cauchy problem for this 
upper triangular system is 



C t = aiE 1 +ba 2 E 2 , 



(2.23) 



C p = (X 2 — Ai) a 2 E 2 , 



with 



<*\ = hb 



bk 



pb 



X 2 — Ai ' 



&2 



kpb 



X 2 — Ai 



and £1, ,£2 defined as in (2.10). 



Finally, we consider the diagonal case b = 0, c = 0, i.e. 



-a 
-d 



Then, it is straightforward to obtain 

C t = hb E u C p = £^2- 

and 



Cu = -kpb E 2 + / C 6 (r) dr. 

./o 



(2.24) 



(2.25) 



(2.26) 



(2.27) 



(2.28) 



3. Solution of the inverse problem. The model equations obtained in the 
previous section describe the time behavior of the tracer concentration in the three 
compartments of the renal system, given the TAC and the transmission coefficients. 
Given such equations, compart mental analysis requires the determination of the tracer 
coefficients by 

1. utilizing measurements of the tracer concentrations provided by nuclear imag- 
ing, and 

2. applying an optimization scheme for the solution of the inverse problem. 
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In nuclear imaging experiments, the reconstructed images can provide information 
on the tracer concentration in the kidneys and in the bladder as well as in the input 
arterial blood as measured in the left ventricle. Specifically, an acquisition sequence 
is set up providing count data sets collected at subsequent time intervals. For each 
data set, an image reconstruction algorithm is applied, ROIs are drawn within the left 
ventricle, the kidneys and the bladder and the corresponding tracer concentrations 
are computed. Obviously, the tracer concentration in the kidneys is an estimate 
of Ct + C p plus the tracer carried by the blood contained in the kidneys' vascular 
system. This last term cannot be identified in the images and will be accounted for in 
the optimization procedure. Specifically, the optimization scheme we will implement 
minimizes, at each time point, the functional 

C := \\(C t + C p ) - C exp \\ 2 + \\C U - C u \\ 2 , (3.1) 

where C n , Ct and C p are the analytical solutions of the direct problem computed at 
the given time point; C u is the concentration measured from the ROI in the bladder 
at the time point; C epx := Ck -[~_ V y^ h •> where is the concentration measured from 
the ROI on the kidneys, is the value of the TAC at the specific time value and V& 
measures the blood fraction with which the kidneys' vascular system is supplied. In 
the following we will assume V& = 0.2, which is a physiologically sound value pQ. 

The minimization of functional C is realized by means of an Ant Colony Optimiza- 
tion (ACO) scheme [3]. ACO is a statistical-based optimization method developed in 
the 1990s with the aim of providing a reliable although not optimal solution to some 
non-deterministic polynomial-time hard combinatorial optimization problems. More 
recently ACO has been generalized to continuos domains and has been applied to a 
wide range of problems. ACO takes inspiration from the way in which ants find and 
carry food to their nest. While an ant is going back to the nest after having taken 
some food, it releases a pheromone trace that serves as a trail for next ants, which 
are able to reach food detecting pheromone. Since the pheromone decays in time, 
its density is higher if the path to food is shorter and more crowded; on the other 
hand, more pheromone attracts more ants and at the end all ants follow the same 
trail. This behavior is paraphrased in ACO identifying the cost functional T with 
the length of the path to food, and the pheromone traces with a probability density 
which is updated at each iteration depending on the value of the cost function for a 
set of states. In practice, at each iteration, the cost function is evaluated on a set of 
P admissible states, and the states are ordered according to increasing values of the 
cost function. Then ACO defines a probability distribution which is more dense in 
correspondence with the cheaper states and, on its basis, Q new states are extracted. 
A comparison procedure identifies the new best P states which form the next set of 
states. Formally the starting point of the algorithm is a set of P states 

B := {U fc = Oi,/c, . . . , u Nik )}, (3.2) 

such that 

V k eScR N , fc = l,...,P (3.3) 

that are ordered in terms of growing cost, namely, .F(Ui) < • • • < J 7 (Up). Next, for 
each j = 1, . . . , N and i = 1, . . . , P, one computes the parameters 

£ p 

m i,j — u j,ii s i,j ~ p Y ^ ^ \ U hP ~ U jA-> (3-4) 

p=l 
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and defines the probability density function 



w iNim itj ,8 itj ](t), (3-5) 
i=l 

with i = 1, . . . , P and £, g real positive parameters to be fixed and il^ = A/ji ?g p] (i). By 
sampling S Q times with the procedure generates Q new states Up+i, . . . ,Up + g, 
enlarging the set B to B = {Ui, . . . , Up + q}. If U/^, . . . , U/c Q are the Q states of B 
of greater cost, the updated B is defined as 

B = B\{V kl ,...,\J kQ }. (3.6) 

This procedure converges to an optimal solution of the problem by exploiting the fact 
that the presence of the weights in the definition of gives emphasis to solutions 
of lower costs since w\ > • • • > wp. This fact, associated with the influence that a 
proper choice of parameters £ and q has on the shape of the Gaussian functions, 
determines the way in which the method tunes the impact of the worse and best 
solutions. The algorithm ends when the difference between any two states of B is less 
than a pre-defined quantity or when the maximum allowable number of iterations is 
reached. The initial set B of trial states is chosen by sampling a uniform probability 
distribution. 

The implementation of ACO for the optimization of the exchange coefficients in 
C is based on the following steps: 

1. The four ACO parameters are fixed as follows. P and Q are chosen as in [4 . 
Specifically, P is a multiple of the number of coefficients to optimize plus one 
and 



Q 



1, (3.T) 



where [•] indicates the floor; q and £ are fixed searching for a trade off between 
the risk of a solution space of a too high complexity and the risk of a too high 
computational demand. 

2. The values of the tracer coefficients are initialized to six random numbers 
picked up in the interval (0, 1) (this choice is based on literature [19]). 

3. The ACO procedure is then run using as C u , C t and C p the solutions in 
the general full matrix case. If, during the ACO iterations, the reconstructed 
values for the tracer coefficients become statistically consistent with values for 
which the associated direct problem is described by a triangular or a diagonal 
matrix, the algorithm automatically switches to utilize the corresponding 
solution in the computation of C. 

In the following section we show how this statistics-based compartmental analysis 
works in the case of synthetic data and real measurements recorded by a micro-PET 
system. 

4. Numerical experiments. Compartmental analysis is a valid approach to 
physiological studies of animal models by means of PET data. An 'Albira' micro- 
PET system produced by Carestream Health is currently operational at the IRCCS 
San Martino 1ST, Genova and experiments with mice are currently performed by using 
different tracers, mainly for applications to oncology. In this section we describe the 
performance of our approach to compartmental analysis in the case of synthetic data 
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Fig. 4.1. Results obtained with the following ACO parameters: P = 13 ; Q = 7, q = 0.015, 
£ = 0.4 for 30 runs of the algorithm. Left panel: red line represents C^, blue line represents Ct, 
black line represents C p ; green line represents C u - Right panel: superimposition of synthetic data 
(white dots) and reconstructed confidence strips of concentrations. 
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ktp 
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hb 


kpb 


ground truth 
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0.02 


0.08 


0.3 


0.3 


Mean 


1.0164 


0.0201 


0.0219 


0.0819 


0.3306 


0.3040 


Std 


0.1060 


0.0094 


0.0075 


0.0049 


0.0341 


0.0199 



Table 4.1 



Simulated values of tracer coefficients; average values and standard deviations over 30 runs of 
ACO. 



simulated by mimicking 'Albira' acquisition for FDG-PET experiments. Then we will 
describe the results of data analysis for two real experiments, performed by using 
FDG. 

4.1. Application to synthetic data. In order to produce the synthetic data 
we started from six initial values for the tracer coefficients. These selected values 



generate a matrix A and a vector K as in (2.4) and (2.5). The corresponding solution 
for Ct, C p and C u associated to the Cauchy problem (2.1) are given in (2.8), (2.9) and 



(2.16). The solutions, represented in Figure 4.1 , left panel, are sampled at 27 time 



points corresponding to the distribution of acquisition times experimentally performed 
by 'Albira'. The red line in the panel represents the TAC Cb that has been obtained 
by fitting with a gamma variate function a set of real measurements acquired from a 
healthy mouse in a very controlled experiment 5|. The vectors corresponding to the 
discretization of C t and C p are then summed together and affected by Poisson noise; 
the same kind of noise is applied to the vector corresponding to the discretization 
of C u . ACO is applied 30 times to these synthetic data set to obtain 30 sets of 
reconstructed values of the exchange coefficients; in correspondence with each set we 
then computed the concentration curves and super-imposed them in Figure |4.1[ right 
panel, thus producing confidence strips for the reconstructed concentrations. In Table 
4.1| we present the average values of the tracer coefficients and the standard deviations 
of the values reconstructed by the 30 runs of ACO over the same set of input vectors. 
The same kind of experiment has been performed again, this time choosing as 



initial values the ones that give rise to an upper triangular system matrix as in (2.22). 



The results for this case are in Figure 4.2 and Table 4.2 The reason why, in this case, 
ACO gives zero to the average and standard deviation of k tp is due to the fact that 
our implementation contains a threshold for the parameter outputs (equal to 10 -3 ). 
In this experiment k tp is under-threshold for all runs and therefore the average and 
standard deviation are set to null. 
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Fig. 4.2. Results obtained with the following ACO parameters: P = 13 ; Q = 7 , q = 0.015, 
£ = 0.4 for 30 runs of the algorithm. Left panel: red line represents C^, blue line represents Ct, 
black line represent C p while green line represent C u • Right panel: superimposition of synthetic data 
(white dots) and reconstructed confidence strips of concentrations. 
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kpb 
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0.4 


0.2 


Mean 


0.8921 





0.0209 


0.1089 


0.4121 


0.1979 


Std 


0.1078 





0.0044 


0.0102 


0.0526 


0.0166 



Table 4.2 



Simulated values of tracer coefficients; average values and standard deviations over 30 runs of 
ACO. 



We agree that in these tests the procedure for generating the synthetic concen- 
trations and the one for reconstructing the tracer coefficients from them are based on 
the same equations (in a sort of 'inverse crime' procedure). However the synthetic 
data are affected by Poisson noise and in any case the aim of these two numerical 
applications was simply to validate the reliability and stability of ACO when applied, 
for the first time, to a compartmental analysis problem. 

4.2. Application to real measurements. For the first experiment we con- 
sidered a healthy murine model injected with FDG and acquired the corresponding 
activity by means of a dynamic acquisition paradigm over 27 experimental time points. 
The images have been reconstructed by applying an Expectation-Maximization iter- 
ative algorithm [2] and ROIs have been drawn on the reconstructed images around 
the left ventricle to reproduce the time activity curve. ROIs have been also drawn 
around the kidneys and the bladder in order to compute the input concentrations. In 



Figure 4.3 , left panel, the red line describes the TAC C5 (we have plotted the solid 
line connecting the measured concentrations in order to distinguish the input function 
from the other concentrations). The green points represent the concentrations for the 
bladder. In order to compute the error bars for the measurements, the concentrations 
have been multiplied times the volume of the bladder and the resulting radioactivity 
is characterized by Poisson statistics, so that an estimate of the standard deviation is 
provided by the square root of the activity. The error bars of the concentration are 
therefore obtained by normalizing this square root with volume. The blue points cor- 
respond to the measured concentration in the kidneys and the error bars are computed 
as before. Then ACO has been applied 20 times against these data. The initializa- 
tion values for the tracer coefficients are the same for all 20 runs and are obtained by 
drawing randomly a 6-ple within (0, l) 6 . The results of this analysis are given in Table 



4.3 containing the average and standard deviation (over the 20 realizations) of the 
tracer coefficients. Then, for each run, the set of reconstructed tracer coefficients are 
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k U p 


hb 


kpb 


initial values 


0.9572 


0.4854 


0.8003 


0.1419 


0.4218 


0.9157 


mean 


1.10 


0.03 


0.04 


0.30 


0.26 


0.25 


std 


0.41 


0.02 


0.03 


0.06 


0.02 


0.02 



Table 4.3 

Initial random values of tracer coefficients; average values and standard deviations over 20 runs 
ofACO. 




Fig. 4.3. Results obtained with the following ACO parameters: P = 25 ; Q = 13 ; q = 0.0001, 
£ = 0.65 for 20 runs of ACO. Left panel: red line represents C^, blue dots Ct + C p while green 
dots represent C u - Right panel: superimposition of concentrations in the bladder (green) and in the 
kidneys (blue) computed by solving the forward problem where the tracer coefficients are reconstructed 
by ACO. 



used to solve the direct problem in order to obtain reconstructions of C u and C p + Ct. 
The confidence strip resulting from the superposition of the 20 reconstructions of the 
concentration are represented in Figure [Oj right panel. A similar analysis has been 
performed in the case of a murine model affected by breast cancer. The results of this 



second case are given in Figure 4.4 and in Table |4~4] 



5. Comments and conclusions. This paper deals with renal flow of a radioac- 
tive tracer, [ 18 F]-FDG, injected into a mouse. The time evolution of tracer concentra- 
tions inside kidney and from kidney to bladder has been modeled by a linear system 
of ordinary differential equations with constant coefficients. The time variation of 
the total concentration of activity inside kidney and bladder (essentially, the sum of 
the solutions) has been estimated through an analysis of micro-PET data. The six 
constant exchange coefficients, which provide information on system's physiology and 
metabolism, have been regarded as unknowns. The related inverse problem has been 
solved by applying a statistical optimization procedure, referred to by the acronym 
ACO. Resulting applications to real and synthetic data have been shown and dis- 
cussed. 

The mathematical approach described in this paper provides estimates of the 
six unknown coefficients. Unlike techniques based on graphical analysis, it does not 
require any distinction between irreversible or reversible uptake of tracer, nor identifi- 
cation of a time value after which suitable expressions evaluated from the data become 
linear in time [18]. Moreover, the graphical methods provide fewer parameters, usu- 
ally slopes and intercepts, which can be interpreted as functions of the original model 
parameters [13]. Technically, the general character of the optimization procedure 
based on ACO makes it applicable to compart mental model structures of a variety 
of types, provided measured data on the related total concentrations of activity are 
available. For example, no practical restriction on the number of compartments in- 
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Fig. 4.4. Results obtained with the following ACO parameters: P = 25 ; Q = 13, q = 0.001 7 
£ = 0.45 for 20 runs of ACO. Left panel: red line represents C^, blue dots Ct + C p while green 
dots represent C u - Right panel: superimposition of concentrations in the bladder (green) and in the 
kidneys (blue) computed by solving the forward problem where the tracer coefficients are reconstructed 
by ACO. 
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0.2503 
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0.19 
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0.05 
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0.01 


0.01 


0.02 



Table 4.4 

Initial random values of tracer coefficients and their average values and standard deviations 
over 20 runs of ACO. 



volved is required if the direct problems involved in the ACO procedure are solved 
numerically, instead of finding the analytic representation of the solutions, as we have 
done in the present paper. Similarly, the fully numeric approach does not require 
limitations such as the constraint that the system matrix is sign-symmetric [8] [17] 
(indeed, the upper and lower triangular matrices considered in Section 2 do not com- 
ply with sign-simmetry) . We have also seen that application of ACO to synthetic data 
corresponding to an upper triangular matrix leads to the final estimate k tp = 0. This 
shows that the ACO approach is also capable of recovering vanishing rate constants. 

The physiological basis for this study relies on the broad utilization of FDG in 
the diagnosing and staging of cancer. In fasting patients, this tracer accurately maps 
the insulin independent glucose metabolism as an index of aggressiveness and growth 
rate of neoplastic lesions. From a theoretical point of view, dynamic PET imaging 
might provide a quantitative estimation of glucose consumption in cancer by mea- 
suring tracer rate constants and serum glucose levels. However, the long acquisition 
time requested by this approach conflicts with the routine activity since it cannot be 
applied to the whole body, it is hardly tolerated by patients and it markedly reduces 
the number of diagnostic procedures. To at least partially overcome this limitation, 
current clinical protocols imply that the whole body is imaged in a single shot late 
after tracer injection. Under this condition, glucose consumption is approximated by 
the measurement of tracer accumulation, normalizing for administered radioactivity 
and distribution mass. However, although this standardized uptake value (SUV) is 
usually considered an adequate surrogate of cancer glucose uptake [21] it also and 
equally depends upon tracer bioavailability. Most often this latter variable is ne- 
glected since blood FDG clearance by tissues is assumed to be relatively stable in the 
different patients. Nevertheless, differently from glucose, a significant amount of FDG 
is excreted in the urines. This tracer sequestration has been already documented and 
attributed to the low affinity of glucose symporters dedicated to glucose reabsorption 
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from preurine in the kidney tubules [TiJ [16]. Our microPET approach offered the 
unique advantage to simultaneously measure the time concentration curves defining 
arterial input function (measured in the left ventricle), kidney tracer handling and 
urinary sequestration in the bladder without any need to evaluate the complex FDG 
transit in renal pelvis [12] and its uncertainties. This direct estimation documented 
that tracer exchange from preurine to kidney parenchima is extremely low. Kidney 
glucose reabsorption is an active process by two Na+-dependent symporters present 
in the apical membrane of tubular cells (SGLT-1 and SGLT-2). The hydroxyl groups 
in D-glucose play a critical role in the recognition and transportation of glucose out 
of the proximal tubules of the kidney. Removal of the hydroxyl group at the 2 posi- 
tion of D-glucose diminishes the affinity of this molecule for both the forms of SGLT 
[TT] . Urinary excretion has been most often considered as a limitation for FDG in 
studying cancer of the urinary tract [9]. However, the evidence of an extremely low 
reabsorption of tracer filtered in the preurine has a more general relevance: it implies 
that tracer loss caused by glomerular filtration markedly affects FDG persistence in 
the blood. Accordingly, SUV might largely differ in cancer lesions characterized by a 
similar metabolism being lower in patients with preserved glomerular filtration rate 
and higher in the presence of an impaired renal function. The disposal of methods 
able to estimate the relationship between the easily available glomerular filtration rate 
and FDG bioavailability might thus improve the clinical standardization of tracer up- 
take as a marker for tumor aggressiveness, beyond the current standard application 
of SUV. 
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